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RESEARCH  BACKGROUND 


The  thermal  stabilization  of  proteins  is  an  intriguing  scientific  problem  with  important 
evolutionary  and  biotechnological  implications  (1-3).  What  mutation(s)  should  be  engineered  to 
increase  the  heat  tolerance  of  a  protein?  There  may  be  no  unique  answer  to  this  question.  Structural 
studies  of  thermophilic  proteins  and  their  mesophilic  counterparts  have  indicated  many  small 
contributions  but  no  single  structural  feature  that  could  account  for  the  increased  thermostability  (4- 
9).  A  powerful  way  to  stabilize  oligomeric  proteins  is  to  increase  the  number  of  intersubunit 
interactions.  However,  in  a  hyperthermophilic  organism  all  proteins  have  to  be  heat  resistant,  even 
monomeric  ones.  Factors  known  to  contribute  to  thermal  stability  include  more  hydrophobic 
contacts  in  solvent-inaccessible  areas,  increased  number  of  salt  bridges,  more  hydrogen  bonds, 
more  proline  residues  and  a  higher  degree  of  compactness. 

In  an  attempt  to  broaden  our  understanding  of  the  principles  of  protein  adaptation  to  high 
temperatures,  we  have  undertaken  the  crystal  structure  determination  of  an  alcohol  dehydrogenase 
(TBAD)  from  Thermoanaerobacter  brockii,  a  bacterium  isolated  from  hot  springs  in  Yellowstone 
National  Park.  TBAD  is  an  NADP  (nicotinamide  adenine  dinucleotide  phosphate)  dependent 
enzyme,  which  reversibly  catalyzes  the  oxidation-reduction  reactions  of  alcohols  to  their 
corresponding  aldehydes  or  ketones  (10).  The  best  substrates  are  secondary  alcohols.  The 
temperature  at  which  half  the  activity  is  lost  after  one  hour  of  incubation  is  93  °C,  and  the  melting 
temperature  is  98°C  (11).  Another  unusual  property  of  this  extremophilic  enzyme  is  its  high 
tolerance  towards  organic  solvents.  TBAD  has  been  used  in  organic  synthesis  as  a  biocatalyst  for 
the  stereospecific  reduction  of  a  broad  range  of  ketones  to  their  corresponding  secondary  alcohols 
(12,  13).  TBAD  is  a  tetrameric  enzyme  of  identical  37,652  Da  subunits  composed  of  352  amino 
acids  (14).  In  contrast  to  eukaryotic  alcohol  dehydrogenases,  which  contain  a  catalytic  Zn2+  as  well 
as  a  structural  Zn2+,  TBAD  has  only  a  single  Zn2+  ion  per  subunit,  reported  to  be  of  catalytic 
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function  (15).  The  ligands  to  the  Zn2+  have  been  identified  by  alanine  mutagenesis  as  Cys37,  His59, 
and  Asp  150  (15). 

The  crystal  structures  of  several  eukaryotic  alcohol  dehydrogenases  have  been  determined, 
including  those  from  cod  and  horse  liver  (CLAD  and  HLAD)  (16, 17),  as  well  as  three  human 
recombinant  forms  (18-20).  The  ligands  to  the  catalytic  Zn2+  in  eukaryotic  alcohol  dehydrogenases 
are  a  histidine  residue  and  two  cysteine  residues.  The  fourth  ligand  is  usually  a  water  molecule  in 
the  absence  of  substrate.  When  the  substrate  is  bound  to  the  enzyme,  the  alcohol  oxygen  or 
carbonyl  oxygen  atom  appears  to  coordinate  directly  to  the  Zn2+,  displacing  the  molecule  of  water. 
Studies  of  the  pH  dependence  of  binding  of  a  variety  of  alcohols  to  HLAD  indicate  that  the  alcohol 
probably  binds  in  the  form  of  alcoholate  anion  (21).  The  negative  charge  on  the  oxygen  atom  makes 
the  alcoholate  a  better  donor  of  a  hydride  ion  for  the  reduction  of  NAD+,  compared  with  the 
undissociated  alcohol.  The  positively  charged  zinc  ion  stabilizes  the  negative  charge  of  the 
alcoholate  in  HLAD. 

The  sequence  of  TBAD  is  remotely  related  to  those  of  eukaryotic  alcohol  dehydrogenases 
with  a  degree  of  sequence  identity  around  30%.  In  contrast,  TBAD  is  highly  homologous  to  a 
recently  crystallized  prokaryotic  alcohol  dehydrogenase  from  Clostridium  beijerinckii  (CBAD) 
with  a  sequence  identity  of  75%  (22).  CBAD  and  TBAD  differ  in  their  thermal  stability  and  in 
their  resistance  to  denaturing  agents.  The  extent  of  resistance  to  denaturation  in  TBAD  depends  on 
the  presence  of  the  cofactor  (22).  The  temperature  at  which  half  the  activity  is  lost  within  one  hour 
is  67  °C  and  93  °C  for  CBAD  and  TBAD,  respectively  (22).  The  increased  thermostability  in 
TBAD  must  be  localized  to  those  25%  of  the  residues,  which  are  different  in  the  two  enzymes. 

We  report  here  the  crystal  structure  of  TBAD  in  a  binary  substrate  complex  with  sec-butanol 
to  a  resolution  of  3.0  A.  Catalysis  does  not  take  place  in  these  crystals  due  to  the  absence  of 
cofactor. 
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EXPERIMENTAL  PROCEDURES 


Crystallization,  X-ray  data  collection  and  processing 

TBAD  was  obtained  from  the  laboratory  of  Dr.  Yigal  Burstein  at  the  Weizmann  Institute  of 
Science.  TBAD  crystals  were  grown  by  the  vapor  diffusion  method  against  a  well  solution 
containing  10%  PEG-8000  (polyethylene  glycol),  5.8%  (v/v)  MPD  (methyl  pentane  diol),  2.5% 
(v/v)  sec-butanol,  50  mM  magnesium  acetate,  0. 1  M  MES  (2[N-Morpholino]ethanesulfonic  acid) 
buffer,  pH  5.8.  The  crystals  belong  to  space  group  P2[2,2i,  with  unit-cell  dimensions  a=80.45  A, 
b=123.08  A,  c=168.03  A.  There  are  four  molecules  of  TBAD  per  asymmetric  unit.  The  volume- 
to-mass  ratio  is  2.86  A3Da‘,  which  is  close  to  the  average  for  protein  crystals  (23).  The  solvent 
content  is  57%. 

Prior  to  data  collection,  the  crystals  were  dialyzed  overnight  against  a  cryoprotectant  solution 
containing  10%  PEG-8000,  20%  (v/v)  MPD,  2.5%  sec-butanol,  50  mM  magnesium  acetate,  0.1  M 
MES  buffer,  pH  5.8.  X-ray  diffraction  data  were  collected  at  CHESS  beamline  A1  using 
synchrotron  radiation  (^=0.908  A)  at  100  °K,  with  a  crystal-to-detector  distance  of  253  mm.  The 
crystals  diffracted  to  2.99  A.  X-ray  diffraction  data  were  processed  using  program  DENZO  and 
reduced  with  SCALEPACK  (24).  Table  1  shows  the  statistics  of  data  collection  and  processing. 

Structure  determination  and  refinement 

The  structure  of  TBAD  was  determined  by  molecular  replacement  with  the  program  AmoRe 
(25)  within  the  CCP4  program  package  (26),  using  data  in  the  resolution  range  from  15  to  4  A. 
Reflections  with  F>ap  were  used  throughout  the  molecular  replacement  and  subsequent  refinement 

calculations.  The  crystal  structure  of  the  CBAD  monomer  (PDB  code  1KEV)  devoid  of  NADP  and 
water  molecules  was  chosen  as  a  starting  model  (27).  The  rotation  search  solutions  were  somewhat 
noisy.  Translation-function  searches  were  subsequently  carried  out  at  the  same  resolution  range  for 
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the  top  35  rotation  solutions.  Four  of  the  solutions,  ranked  1,  3,  13,  and  26  in  the  rotation  search, 
were  found  to  be  outstanding  in  the  translation  search  with  correlation  coefficients  at  least  38% 
higher  than  the  rest  of  solutions.  The  correctness  of  these  four  solutions  was  validated  by  the  fact 
that  they  pack  well  to  form  the  tetramer,  and  they  conform  to  the  subunit  relationships  revealed  by  a 
self-rotation  search  (data  not  shown).  The  tetramer  has  non-crystallographic  222  symmetry.  The 
model  of  the  tetramer  generated  from  the  four  top  solutions  yielded  an  R-factor  of  36.1%  and  a 
correlation  coefficient  of  63.6%. 

The  molecular-replacement  model  was  subsequently  modified  by  mutating  non-conserved 
residues  to  Ala  and  setting  all  B-factors  to  20  A2.  Refinement  was  carried  out  using  X-PLOR, 
versions  3.1  and  3.8  (28).  Rigid-body  refinement  was  carried  out  to  improve  the  positional  and 
orientational  parameters  of  the  four  subunits.  Strict  non-crystallographic  symmetry  was  maintained 
in  this  and  all  subsequent  refinement  cycles.  A  series  of  twenty-residue  omit  (F0-Fc )  maps  were 
calculated  and  used  to  retrace  the  whole  model  using  the  TBAD  amino  acid  sequence.  The  Zn2+  ion 
and  the  substrate  were  located  from  a  residual  (F0-Fc )  map  after  all  protein  residues  were  accounted 
for.  Manual  rebuilding  between  refinement  cycles  was  performed  with  the  computer-graphics 
software  O  (29).  The  model  was  subjected  to  several  cycles  of  simulated  annealing  by  “heating” 
the  system  to  3500  °K  and  slowly  “cooling”  it  to  300  °K  in  time  steps  of  0.00025  ps.  This  reduced 
the  R-value  to  26.3%.  Positional  parameters  and  temperature-factors  were  then  refined  in 
alternative  cycles.  A  (2F0-FC )  Fourier  map  was  computed  at  every  stage  of  the  refinement  and  the 
model  was  adjusted  manually  to  fit  this  electron-density  map.  Water  molecules  were  not  added, 
except  for  the  one  in  close  proximity  to  the  zinc  ion.  Side  chain  density  was  good  for  all  residues 
except  for  a  small  portion  of  side  chains  on  the  surface.  The  Ramachandran  plot  shows  that  84.6% 
of  the  residues  are  in  most  favored  regions  with  none  in  the  disallowed  regions,  as  calculated  with 
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program  PROCHECK  (30).  The  statistics  of  the  crystallographic  analysis  and  the  final  model  are 
listed  in  Table  2.  The  coordinates  will  be  deposited  at  the  Protein  Data  Bank. 

SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

Structure  description 

As  in  all  other  alcohol  dehydrogenases  of  known  structure,  the  monomer  is  made  up  of  two 
domains,  the  catalytic  domain  and  the  cofactor-binding  domain  (Fig.  1).  The  catalytic  domain 
contains  the  N-  and  C-terminal  segments,  consisting  of  residues  1  to  149  and  297  to  352.  There  are 
two  P-sheets  and  five  a-helices  in  the  catalytic  domain.  Strands  1  and  2  form  a  small  antiparallel  p- 
meander  on  the  surface  of  this  domain,  while  strands  3,4,5,12  and  13  form  a  mixed  P-sheet.  The 
loop  between  strands  3  and  4  is  very  long  and  contains  a  helical  stretch  al.  Following  P5  there  are 
four  a-helical  segments  in  a  row.  The  last  of  these,  a5,  links  the  catalytic  domain  to  the  nucleotide¬ 
binding  domain. 

The  nucleotide-binding  domain  consists  of  the  contiguous  stretch  of  residues  157  to  292.  It 
forms  a  typical  Rossmann  fold  with  six  parallel  P-strands,  P6  to  pi  1.  This  p-sheet  is  flanked  by  a- 
helical  segments,  a6  and  a7  on  one  side  of  the  sheet,  and  a8  and  a9  on  the  other  side  of  the  sheet. 
Following  the  last  strand  in  this  sheet,  pi  1 ,  there  is  a  long  loop  which  crosses  over  back  into  the 
catalytic  domain  to  lead  into  helix  alO,  strand  pi 2,  helix  al  1  and  finally  strand  pi 3.  The  two 
domains  are  separated  by  a  deep  active  site  cleft  that  is  accessible  to  the  solvent. 

The  active  site 

The  location  of  sec -butanol  and  the  zinc  ion  in  the  active  site  cleft  was  determined  from  a 
residual  (F0-Fc)  difference  Fourier  map  (Fig.  2).  The  binding  site  for  the  substrate  sec-butanol 
consists  of  a  hydrophobic  region  of  residues  Ile86,  Leul07,  Tyr267,  and  Leu294,  as  well  as  the 
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hydrophilic  residues  His59  and  Asp  150  (Fig.  3).  The  hydroxyl  group  of  the  substrate  is  hydrogen- 
bonded  to  Ne2  of  His59  and  one  of  the  carboxylate  oxygen  atoms  of  Asp  150  with  distances  of  3.2 
and  3.0  A,  respectively.  These  residues  have  been  shown  to  be  essential  for  catalysis  (15). 
Surprisingly,  the  zinc  ion  is  not  in  direct  proximity  to  the  substrate  but  it  is  located  some  9  A  away 
from  the  hydroxyl  oxygen  atom  of  the  substrate.  Residues  surrounding  the  zinc  ion  include  Cys37, 
Thr38,  Met337,  Metl51  and  a  water  molecule  (Table  3).  His59  and  Asp  150  are  not  in  direct  contact 
with  the  zinc  ion.  However,  if  the  zinc  ion  were  to  be  involved  in  catalysis,  a  conformational 
change  would  be  required,  presumably  upon  NADP  binding,  which  would  bring  the  zinc  ion  closer 
to  the  substrate.  Indeed,  substantial  conformational  changes  between  the  apo  and  holo  forms  of  liver 
alcohol  dehydrogenases  have  been  observed  (31). 

Structural  features  that  account  for  the  thermostability 

The  degree  of  sequence  identity  between  TBAD  and  CBAD  is  75%.  Factors  that  account  for 
the  increased  thermostability  in  TBAD  must  therefore  be  located  in  the  25%  of  the  residues  that  are 
different  between  the  two  enzymes  (Fig.  4).  The  sequence  alignment  is  unambiguous  since  there 
are  no  gaps.  There  are  a  total  of  87  residues  that  are  different  between  the  two  enzymes.  As  shown 
in  Table  4,  the  most  significant  differences  in  amino  acid  composition  between  the  two  enzymes  are 
for  alanine,  proline  and  serine. 

Proline  residues — The  number  of  proline  residues  increases  from  13  in  the  mesophilic  enzyme 
to  21  in  the  thermophilic  enzyme  (Table  4).  Proline  reduces  the  main-chain  flexibility  of  the 
backbone,  thereby  decreasing  the  entropy  and  increasing  the  free  energy  of  the  unfolded  state  (32). 
Consequently  any  replacement  to  a  proline  should  destabilize  the  unfolded  state,  thereby  stabilizing 
the  folded  state.  Interestingly,  these  replacements  do  not  occur  at  random  positions  but  at 
preferential  sites.  Of  the  eight  extra  prolines  in  TBAD,  Pro  177,  Pro222  and  Pro316  are  found  at  the 
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N-terminal  position  of  helices.  A  residue  in  this  position  is  partly  outside  the  helix  and  partly 
inside,  i.e.  the  proline  ring  is  not  part  of  the  helix  but  the  carbonyl  oxygen  forms  a  a-helical 
hydrogen  bond  with  the  amide  nitrogen  four  residues  upstream.  Proline  residues  rarely  occur  in 
helix  interiors,  where  they  would  cause  a  kink  in  the  helix  (33,  34).  Side  chain  hydrogen  bonding  to 
the  backbone  in  a  helix  is  rarely  to  the  amide  nitrogen  of  the  first  residue  but  frequently  to  the 
second  or  third  residue.  The  lack  of  a  hydrogen  atom  on  the  amide  nitrogen  of  a  proline  residue  in 
position  one  of  a  helix  seems  therefore  not  be  a  hindrance.  Single  proline  site  directed  mutagenesis 
of  CBAD  revealed  the  following  changes  in  the  melting  temperature:  an  increase  for  the  alanine-to- 
proline  substitution  at  position  177  and  the  leucine-to-proline  substitution  at  position  316  of  0.5  °C 
and  10.8  °C,  respectively;  and  a  decrease  for  the  histidine-to-proline  substitution  at  position  222  of 
6.5  °C  (1 1).  Position  222  is  on  the  surface  of  the  protein.  Replacing  a  histidine  for  a  proline  at  this 
position  leads  to  a  loss  of  solvation  energy,  which  would  account  for  the  observed  destabilization. 
Pro316  is  the  first  residue  in  a  4-residue  long  stretch  of  310  helical  conformation,  linking  an  a-helix 
to  a  p-strand.  Replacing  a  bulky  leucine  for  the  smaller  proline  may  relieve  some  overcrowding  in 
this  region  containing  Metl47,  Metl52  and  Val309.  Therefore,  one  may  conclude  that  the 
introduction  of  a  proline  at  the  beginning  of  a  helix  does  not  stabilize  the  protein  unless  it  is  a 
special  position.  Pro24  is  located  at  position  two  in  a  P-tum  (Fig.  5),  a  position  that  is  frequently 
populated  by  prolines.  The  serine-to-proline  mutant  of  CBAD  at  this  position  has  a  3.9°C  higher 
melting  temperature  (11).  The  introduction  of  a  kink  at  this  position  is  facilitated  by  the  fixed 
angle  in  proline. 

Alanine  and  serine  residues — As  shown  in  Table  4,  TBAD  has  eight  more  Ala  residues 
compared  to  CBAD.  Three  of  the  replacements  to  alanine  occur  in  a-helical  regions.  An  alanine 
increases  the  stability  of  a  helix  because  of  the  unrestricted  rotational  freedom  of  the  methyl  side 
chain  (35,  36).  Introduction  of  alanine  residues  at  other  positions  involve  replacement  of  polar  side 
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chains  such  as  serine  and  arginine  in  buried  regions  (Ser35,  Serl68,  Ser246,  Ser250,  Ser263, 
Arg238),  as  well  as  replacement  of  valine  with  the  smaller  and  less  hydrophobic  alanine  on  the 
surface  of  the  protein  at  position  32  (Fig.  4).  These  replacements  also  provide  entropic 
stabilization.  Some  serine  residues  in  the  core  of  the  protein  are  replaced  with  the  more  apolar 
threonine  residues  without  loss  of  hydrogen  bond  capacity  (Ser92,  Serl54,  and  Serl69).  Two 
exposed  serine  residues  at  positions  254  and  270  are  replaced  with  the  more  hydrophilic  lysine  and 
glutamic  acid,  respectively. 

Surface  characteristics  and  subunit  interfaces — TBAD  is  more  compact  than  its  mesophilic 
counterpart.  The  accessible  surface  area  in  TBAD  is  7.7%  smaller  than  in  CBAD  (Table  5).  Tighter 
packing  has  indeed  been  suggested  as  a  mechanism  for  enhanced  thermostability  in  proteins  (4). 
Moreover,  the  area  buried  upon  oligomerization  is  significantly  larger  in  TBAD  (Table  5);  i.e.  there 
are  more  extensive  contacts  between  the  subunits  in  the  thermophilic  enzyme.  These  intermolecular 
contacts  include  salt  bridges,  hydrogen  bonds  and  hydrophobic  interactions. 

Salt  bridges,  hydrogen  bonds  and  van  der  Waals  contacts — A  major  cause  of  higher 
thermostability  can  be  attributed  to  an  increase  in  the  number  of  van  der  Waals  contacts,  hydrogen 
bonds  and  salt  bridges,  particularly  those  at  subunit  interfaces  (Table  6).  Although  the  TBAD 
monomer  has  fewer  hydrogen  bonds  in  comparison  to  CBAD,  due  to  the  loss  of  buried  polar 
residues,  this  is  more  than  offset  by  the  creation  of  extra  salt  bridges  and  van  der  Waals  contacts 
both  within  and  between  the  monomers.  Furthermore,  additional  intermolecular  hydrogen  bonds  are 
gained  in  TBAD.  Salt  bridges  are  important  stabilizing  forces  in  proteins  if  they  occur  in  regions 
inaccessible  to  the  solvent  (37-39).  The  strength  of  a  salt  bridge  depends  on  the  environment,  i.e. 
the  local  dielectric  constant,  and  the  bond  distance.  There  are  3 1  salt  bridges  in  TBAD  versus  24  in 
CBAD  (Table  6).  Some  of  the  salt  bridges  in  TBAD  are  involved  in  extensive  networks.  The 
energy  gain  in  a  network  is  larger  than  the  sum  of  the  pairwise  interactions  between  the  constituent 
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charges  (40).  An  example  of  a  sequestered  network  of  six  salt  bridges  -  some  intramolecular  and 
some  intermolecular  -  is  depicted  in  Fig.  6.  Only  two  of  these  salt  bridges  exist  in  CBAD. 
Hydrogen  bonds  occur  between  a  hydrogen  bond  donor  and  acceptor  pair  with  a  distance  of  less 
than  3.3  A  at  an  angle  greater  than  90  °(4).  There  are  218  intermolecular  hydrogen  bonds  in  TBAD 
versus  167  in  CBAD  (Table  6).  Accurate  accounting  of  salt  bridges  and  hydrogen  bonds  in  TBAD 
is  difficult,  as  the  crystal  structure  was  determined  at  3.0  A.  A  few  genuine  interactions  may  have 
been  missed,  and  the  validity  of  other  assigned  bonds  may  be  questioned  due  to  large  interatomic 
distances.  However,  significant  differences  in  the  number  of  interactions  between  the  thermophilic 
and  the  mesophilic  enzyme  clearly  suggest  mechanisms  to  enhance  the  thermal  stability  of  TBAD. 
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TABLE  I.  Statistics  of  Data  Collection  and  Processing 


Resolution 

(A) 

Reflections 

Completeness 

all  3.10-2.99 

redundancy 

all  3.10-2.99 

<I/<T> 

all  3.10-2.99 

Emerge 

2.99 

30323 

0.863  0.890 

4.4  4.8 

12.1  4.6 

0.083 

*Rmerge  =Zh£ilIih"<Ih>l®h^<Ih>  where  <Ih>  is  the  mean  intensity  of  the  i  observations  of  reflection  h. 

TABLE  II  Refinement  Statistics 


Resolution  (A) 

50-2.99 

Reflections  (work/free) 

26756/2986 

R*/Rfree** 

0.211/0.265 

No.  of  protein  atoms 

2638 

r.m.s.  bond  lengths  (A) 

0.02 

r.m.s.  bond  angles  (°) 

1.89 

Avg.  B-factors  (main/side)  (A2) 

19.53/23.15 

£IF0-FCI 

*  R=  Fo,  observed  structure  factor;  Fc  calculated 

£F0 

structure  factor 

**  Rfree,  defined  as  in  R  for  a  randomly  selected  set  of  10%  of  the 
reflections 
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TABLE  III  Distances  between  the  zinc  ion  and  surrounding  residues  (A) 


Cys37  Sy 

3.2 

Thr38  Oyl 

4.4 

Met337  S5 

3.6 

Water  0 

3.0 

Metl5 1  S5 

4.9 

TABLE  IV  Differences  in  the  amino  acid  composition 


TBAD 

CBAD 

TBAD -CBAD 

Ala 

35  (.099) 

27  (.077) 

8  (0.022) 

Arg 

14  (.040) 

14  (.040) 

0 

Asn 

10  (.028) 

13  (.037) 

-3  (-0.009) 

Asp 

20  (.057) 

21  (0.060) 

-1  (-0.003) 

Cys 

4  (0.011) 

5  (0.014) 

-1  (-0.003) 

Gin 

3  (0.009) 

6  (0.017) 

-3  (-0.008) 

Glu 

21  (0.060) 

18  (0.051) 

3  (0.009) 

Gly 

43  (0.122) 

45  (0.128) 

-2  (-0.006) 

His 

10  (0.028) 

11  (0.03.1) 

-1  (-0.003) 

He 

26  (0.074) 

26  (0.074) 

0 

Leu 

23  (0.065) 

28  (0.080) 

-5  (-0.015) 

Lys 

24  (0.068) 

21  (0.060) 

3  (0.008) 

Met 

15  (0.043) 

18  (0.051) 

-3  (-0.008) 

Phe 

14  (0.040) 

11  (0.031) 

3  (0.009) 

Pro 

21  (0.060) 

13  (0.037) 

8  (0.023) 

Ser 

9  (0.026) 

17  (0.048) 

-8  (-0.022) 

Thr 

13  (0.037) 

11  (0.031) 

2  (0.008) 

Trp 

4  (0.011) 

4  (0.011) 

0 

Tyr 

6  (0.017) 

7  (0.020) 

-1  (-0.003) 

Val 

37  (0.105) 

35  (0.100) 

2  (0.005) 

The  amino  acid  composition  is  based  on  SWISSPROT  entries  P14941  and  PP25984  for  TBAD  and  CBAD,  respectively.  Values  in 
parentheses  represent  mol  fractions.  Differences  larger  than  one  percent  are  shown  in  bold. 
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TABLE  V  Surface  area  of  the  tetramer 


TBAD 

CBAD 

Difference  (%) 

Total  accessible  surface  area  (A2) 

44996 

48775 

-7.7 

Area  buried  in  tetramer  interfaces  (A2) 

14684 

14397 

2.0 

Fraction  of  buried  area 

0.326 

0.295 

-3.1 

Percentage  of  exposed  hydrophilic  residues 

71.4 

69.2 

2.2 

Percentage  of  exposed  hydrophobic  residues 

28.6 

30.8 

-2.2 

Percentage  of  buried  hydrophilic  residues 

19.4 

24.5 

-5.1 

Percentage  of  buried  hydrophobic  residues 

80.6 

75.5 

5.1 

TABLE  VI  Interactions  in  CBAD  and  TB  AD 


Intermolecular  contacts 

Intramolecular  contacts 

H-bonds 

salt  bridges 

van  der  Waals 

H-bonds 

salt  bridges  van  der  Waals 

CBAD 

167 

8 

774 

538 

16 

4383 

TBAD 

218 

13 

862 

476 

18 

4689 

Difference 

51 

5 

88 

-62 

2 

306 

The  contacts  were  calculated  with  program  CONTACT  within  CCP4  crystallographic  suite  of  programs  (26).  Threshold 
values:  salt  bridges  4.0,  H-bonds  3.4,  van  der  Waals  contacts  4.0  A. 
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FIG.  1.  Ribbon  drawing  of  the  TBAD  monomer.  Helices  and  strands  are  shown  in  dark  and  light 
gray,  respectively.  The  enzyme  is  made  up  of  two  domains,  the  nucleotide-binding 
domain  and  the  catalytic  domain.  The  two  domains  are  separated  by  a  deep  active  site  cleft 
that  is  freely  accessible  to  the  solvent.  The  Zn2+  ion  is  located  at  the  bottom  of  the  cleft  in 
the  catalytic  domain.  This  figure  was  created  with  programs  MOLSCRIPT  (41)  and 
RASTER3D  (42). 

FIG.  2.  Residual  (F0  -  Fc)  difference  Fourier  map  contoured  at  a  level  of  four  standard  deviations 
above  the  mean.  Fc  contains  contributions  from  all  protein  atoms  but  no  other  atoms. 

Peakl  has  been  interpreted  as  the  zinc  ion  due  to  its  proximity  to  Cys37.  Peak2  has  been 
interpreted  as  sec-butanol  due  to  its  proximity  to  the  catalytic  residues  His59  and  Asp  150 
and  a  hydrophobic  region  of  residues  Ile86,  Leu  107,  and  Leu294.  These  are  the  two  most 
prominent  peaks  in  this  map.  Residues  in  the  vicinity  of  the  active  site  are  depicted  in  ball- 
and-stick.  This  drawing  was  made  with  program  O  (29). 

FIG.  3.  Ball-and-stick  representation  of  residues  in  the  vicinity  of  the  active  site.  The  binding  site 
for  the  substrate  sec-butanol  consists  of  a  hydrophobic  region  including  residues  Ile86, 
Leu  107,  Tyr267  (not  shown  for  clarity)  and  Leu294,  as  well  as  the  catalytic  residues 
His59  and  Asp  150.  The  zinc  ion  is  not  in  direct  proximity  to  the  substrate  but  it  is  located 
some  9  A  away  from  the  hydroxyl  oxygen  atom  of  the  substrate.  Residues  surrounding  the 
zinc  ion  include  Cys37,  Thr38,  Metl51,  Met337,  and  a  water  molecule. 
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FIG.  4.  Sequence  Alignment  of  TBAD  with  CBAD.  The  lines  marked  “TBAD”  contain  only 
changes  from  CBAD.  Residues  marked  “b”  are  buried  (accessibility  <  25  A2),  residues 
marked  “e”  are  exposed  in  the  tetramer  (accessibility  >  25  A2). 

FIG.  5.  Ball-and-stick  representation  of  a  region  containing  three  proline  residues  in  the 

thermophilic  enzyme  at  positions  20, 22  and  24.  TBAD  is  shown  in  red  and  CBAD  in 
cyan.  Pro20  is  conserved  in  both  enzymes  whereas  residues  22  and  24  in  CBAD  are 
alanine  and  serine,  respectively.  Introduction  of  a  proline  in  position  22  increases  the 
number  of  van  der  Waals  contacts  with  nearby  Ile28.  Residue  24  is  located  at  position  two 
of  a  P-tum,  a  position  that  is  frequently  populated  by  prolines.  There  is  little  or  no  entropy 
loss  during  folding  for  a  proline  at  this  position. 

FIG.  6.  Ball-and-stick  representation  of  the  extensive  network  of  salt  bridges  involving  three 

subunits  in  TBAD.  Residues  are  labeled  red,  green,  and  white,  according  to  the  subunit 
they  belong  to.  R97  from  the  “white”  subunit  forms  an  intramolecular  salt  bridge  with  E94 
(shortest  distance  3.4A).  Across  the  subunit  interface  from  R97  the  carboxylate  of  D237 
from  the  “red”  subunit  is  located  at  a  distance  of  5.8  A,  normally  considered  too  long  for  a 
salt  bridge.  However,  in  a  milieu  of  low  dielectric  constant  these  two  opposite  charges 
might  still  exert  some  attractive  force  on  each  other.  This  same  residue  D237  also  interacts 
with  K257  from  the  same  subunit  as  well  as  with  R304  from  the  “green”  subunit  at 
distances  of  3.1  and  3.4  A,  respectively.  R304  from  the  “green”  subunit  in  turn  interacts 
with  El 65  from  the  “red”  subunit  at  a  distance  of  4.1  A.  The  carboxylate  of  El  65  from  the 
“red”  subunit  is  7.1  A  from  the  guanido  group  of  R301  from  the  “green”  subunit,  too  far  to 
be  considered  a  salt  bridge.  Finally,  this  R301  interacts  with  El 60  from  the  same  “green” 
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subunit  at  a  shortest  distance  of  2.5  A.  There  are  five  salt  bridges  depicted  in  this  figure  - 
six  if  you  count  the  long  intermolecular  bridge  between  R97  and  D237.  Of  these  only  two 
are  present  in  mesophilic  CBAD. 


19 


1  MKGFAMLGINKLGWIEKERPVAGSYDAIVRPLAVSPCTSDIHTVFEGALGDRKNMILGHE 


ffi  0 
ft  0) 
H  0 


ft  (D 
ft  0 

ft  0 
C  d) 

ft  0) 


>  45 
U  0 
W  0 


99 

pq  pq 

U  Eh 


Eq 

0 

&  ft  0 

c 

H 

>  43 

> 

W 

<3  45 

<3 

o 

£ 

rH  rH 

rH  rH 

^ — I 

rH 

VO  VO 

CN  CN 

00 

00 

rH  ^ — 1 

X — \ 

rH 

99 

PQ  PQ 
U  Eh 


pq  .  . 

U  Eh 


99 

pq  pq 

U  Eh 


J  0 


CM  CM 


pq 

O  Eh 


i  t=U  0 
h3 
H 


ft  43 


H  45 


S0 
43 

ft  d) 
£  0 


ft  d) 
ft  0 


w 

1-3  ft  43 


ft  d) 
ft  0 
J  43 
H  43 
ft  0 
£  43 


r- I  rH 
O  O 
00  00 


99 

se 


viT 


V21V 


r 


6 


LIST  OF  ALL  PUBLICATIONS  AND  TECHNICAL  REPORTS 

Zhang,  Z.,  Djebli,  A.,  Shoham,  M.,  Frolow,  F.,  Peretz,  M.,  Burstein  Y.,  Crystal 
parameters  of  an  alcohol  dehydrogenase  from  the  extreme  thermophile 
Thermoanaerobium  brockii.  J.  Mol.  Biol.  230,353-355  (1993). 

Shoham,  M.,  Protein  Adaptation  to  Extreme  Envionments.  Army  Research  Office 
Biosciences  Workshop,  Cashiers,  North  Carolina,  12-14  May  (1997). 

Shoham,  M.,  Protein  Adaptation  to  Extreme  Environments.  Yellowstone  National  Park 
125th  Anniversary  Symposium  on  Life  in  Extreme  Environments.  Montana  State 
University,  Bozeman,  Montana,  16-18  May  (1998). 

LI,  C.,  Sandriyana,  S.,  H.,  Joel,  and  M.,  Shoham,  Crystal  structure  of  a  thermophilic 
alcohol  dehydrogenase  substrate  complex  from  Thermoanaerobacter  brockii  and 
comparison  to  a  mesophilic  counterpart  Submitted  for  publication 


26 


LIST  OF  ALL  PARTICIPATING  SCIENTIFIC  PERSONNELS 


Menachem  Shoham,  P.  I.,  Associate  Professor  of  Biochemistry 

Chunmin  Li,  Postdoctoral  fellow 
Sandriyana  Soelaiman,  M.  S.  student 
Joel  Heatwole,  undergraduate  student 


27 


